Fractional Vector Calculus and Fractional Special Function 
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Fractional vector calculus is discussed in the spherical coordinate framework. A variation of the 
Legendre equation and fractional Bessel equation are solved by series expansion and numerically. 
Finally, we generalize the hypergeometric functions. 

PACS numbers: 02.30.Gp, 45.10.Hj, 02.40.Ky 



O 

(N 
C 

1—5 



43 

Oh- 

43 

•*-> 
c3 



> 
00 

oo 
(N 

o 

O 



I. INTRODUCTION 

Fractional calculus is the calculus of differentiation and 
integration of non- integer orders [3, 0] ■ During last three 
decades, fractional calculus has gained much attention 
due to its demonstrated applications in various fields of 
science and engineering, such as anomalous diffusion 
fractional dynamical systems fractional quantum 

mechanics fractional statistics and thermodynamics 
Q, to name a few. 

Fractional vector calculus (FVC) is important in de- 
scribing processes in fractal media, fractional electrody- 
namics and fractional hydrodynamics [H, fXoj | . But an 
effective FVC is still lacking. There are many problems 
in defining an effective FVC. One is that fractional in- 
tegral and fractional derivative are defined "half , that 
is to say, they are defined only on the right or the left 
side of an initial point. And if we make fractional series 
expansion, functions are all expanded on the right or the 
left neighborhood of the initial point. We cannot across 
this cutting point. So if we want to describe the behavior 
near the initial point, we need define both the right and 
the left functions. The situation will become even more 
complicated when we deal with high dimensions. In this 
letter we will define FVC in spherical coordinate frame- 
work. Since in the spherical coordinate framework the 
radius is naturally bounded to the positive half, we need 
just one fractional derivative. 

Using this frame, we will discuss the Laplacian equa- 
tion with fractional radius derivative in 3-d space and 
the heat conduct equation in 2-d space with cylin- 
drical symmetry. As a result, the corresponding special 
functions will be generalized. Finally, we will generalize 
hypergeometric functions fl2l - [l4j |. 



II. FRACTIONAL CALCULUS 

To start, let's briefly recall some basic facts in frac- 
tional calculus 0, [H. There are many ways to define 
fractional integral and fractional derivative. In this let- 
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ter we will use Riemann-Liouville fractional integral and 
Caputo fractional derivative. 

Let f(x) be a function defined on the right side of a. 
Let a be a positive real. The Riemann-Liouville frac- 
tional integral is defined by 
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r(a) J a 



(s-o a "7(Ode- (!) 



The integral has a memory kernel. 

Let A = [a] + 1. The Caputo fractional derivative is 
defined by 
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They have the following properties: 
D a (x - af = 0, p G {0,1,..., [a]}; (3) 

D a {x - af = T jy + V) A x - a)^" Q ,/3 otherwise; (4) 
v ; r(/3 - a + 1) v ; w 

These properties are just fractional generalizations of 
d n m ml 



dx n ' 



(m — n)\ 



neN, m ^ 0, (6) 



dx x rn = - 

(m+l)\ 



(7) 



For a 'good' function, one can define its fractional Tay- 
lor series 



f(x) = J2(D a ) m f(x)\ x=a -[(n m -l]- (8) 
Explicitly, 



m=0 



(J«) m . 1 = 



1 



T(ma + 1) 



(x-a) ma . 



(9) 
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FIG. 1: Spherical framework with fractional radius deriva- 
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III. FRACTIONAL VECTOR CALCULUS 

As has been aforementioned, fractional derivative is 
defined only on the right or the left half of the real line, 
which gives complications in defining an effective FVC 
with Cartesian coordinates. So it may be more feasible 
to do FVC with spherical coordinates. 

I 



A. Spherical coordinates 

The spherical coordinates of 3-dimension space is a 
triplet (r,9,<j)). In one-order calculus, the gradient of 
a scale function u(r,9, 4>) is 



Grad u 



d 



1 d 



or r o9 



_1 d_ 

' rsinf? d(j) 



(10) 



We generalize this definition to fractional calculus as 



Grad" u = e r D?+e e 



r(a + l) d r(a + l) d 
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A-a-l d 



r a sm a 6 dcj> 

(11) 

u{p, e, (j>)dp. 



where D% u — r ^_ a) fi(r - p) 
For anisotropic space, a is a function of and <fi- F° r 
isotropic space, a is a constant and independent of 9 and 
<j>. We will just consider the isotropic case. 

By this definition, the real space metric is changed to 
an effective metric. This can be easily seen from Fig. 
Q] The radius \OA\ now is [/« • 1]|^ = p^r^; 

the arc length NA — T (a+i) r A@A and the arc length 

LA' = r(a+i) r A sm "^ ( t )A - This kind of metric is not 
addictive since \OB\ ^ \OA\ + \AB\ even when O, A 
and B arc on the same straight line. This is due to the 
non-locality of the fractional operations. 



By the above generalization, the divergence of a vector function A = (A r , Ag, A^) is 

div " A = 2a 1 an D?(r 2a sm a dA r ) + T £ + ? -?-(r a S m a 0A B ) + ^ £(r"^ ). 
H a sin r Za sm 9 o9 r Za sm a 1 



(12) 



The Laplacian of a scale function u(r, 9, <p) will be 
A a u = div"Grad Q u 



1 



r 2 (a + i) d . . ao d . r 2 (a + i) d 2 



-D ? (r**D? u) + -^-^^(sin"c9- „) + 



sin Q 6» d9 y 09 ' r 2a sin 2a 9 d6 2 



(13) 



The rotor operator can be defined as well. Since in The divergence is 
this letter we do not deal with the rotor, its definition 
will not be given. 



B. Polor coordinates 

Likewise, the fractional gradient operator of two di- 
mensions in polor coordinates can be defined as 



div« A = + V± ^ A §- e Ae. (15) 



The Laplacian of a scale function u(r, 9) is 

A" u = div a Grad Q u 
1 

86 ) ' ~° 



Gr^u=[e r D ? + e 8 T J^ll)u. (14) = ^(^>) + ^^|U (16) 
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IV. FRACTIONAL SPHERICAL EQUATION 
AND FRACTIONAL CYLINDRICAL EQUATION 



A. Fractional Laplacian equation 



In this section, we consider the Laplacian equation 
with the 3-d Laplacian operaror defined above and the 
fractional heat conduct equation in a 2-d space with 
cylindrical symmetry. 



With the Laplacian operator defined above, the Lapla- 
cian equation becomes 



J 



r 2 (a + i) d 



r 2a sm a 9 09 

r 



d_ 

dd 



T 2 (a + 1) d 2 



r 2a sm 2a t 



(17) 



This equation can be solved by separation of variables. 
Let u — R(r)Q(6)$>((j)) and substitute, the result is 



1 d f , a d& 
0sin a dO \ dd 



1 d 2 $ 
$sin 2Q d(b 2 



R 



D?{r 2a D«R) = Ar 2 (a + 1) 



-A, (18) 



(19) 



The second equation can be solved by fractional series 
expansion. Let R = X)m=-oo c raf ma , and substitute, the 
result is that the nonzero components are the ones with 
m satisfying 



T(ma + a + 1) 2 

_ _ = Ar (a + 1). 

1 [ma — a + 1 j 



(20) 



Otherwise, c m =0. 

The first equation is a variation of the ordinary spher- 
ical harmonic equation. By further decomposition, it 
turns to 



d 2 $ 



m 2 $ = 0, (m = 0,l,2,3,...) (21) 



Legendre function 
even, lamda=12 




FIG. 2: Legendre even function. Since the fractional index 
a occurs as a multiplicative factor in the variation version of 
the Legendre equation, small difference from 1 will not make 
large changes to the profile. We calculated the function with 
other small difference a's, but the results are not shown for 
their curves are very close to each other. Notice the direct of 
the curves. 



sin a d ( . an de 



m 2 - Xsm 2a 9. (22) 



The first equation is simple. The second equation 
can be transformed by changing variables x = cosO and 
p(x) = 0(0) to 



(1 



A- 



(1 



r 2\a 



p = 0. (23) 



This is a variation of the ordinary associated Legendre 
equation [0, HI]- By setting a = 1, we recover the ordi- 
nary associated Legendre equation. 

When m — it is the Legendre equation. Make the 
ansatz p(x) = J^^Lo c n x ™ an d substitute, we find that 
the terms with even powers of x and the terms with odd 
powers of x are independent, so we can write p(x) = 



(x) + ci 
coefficients is 



Podd(x). The relation of successive 



-A 



Cn+2 



(n + 2)(n + l)' 



(24) 



Let a = 1 and A = 1(1 + 1), we will get the Legendre poly- 
nomials. We calculated numerically the functions with 
some different values of the parameters and displayed 
the results in Fig. [5] and Fig. [3] 



B. Fractional cylindrical equation 

The heat conduct equation in 2-d space is 
du 



dt 



a A2U. 



(25) 
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FIG. 3: Legendre odd function. Notice the direct of the 
curves. 

Here A 2 is the two dimensional Laplacian operator. Us- 
ing the fractional cylindrical Laplacian operator defined 
above, this equation becomes 

^> ?(r ^„) + If<£^„. (26) 



Of 



de 2 



By separation of variables u(t, r, 9) — R(r)<d(6)T(t), it 
can be decomposed to 



T' + a 2 k 2 T = 0, 



(27) 
(28) 




FIG. 4: Parameters of fractional Bessel function, v as a func- 
tion of a and p. 



Fractional Bessel functions 
v=3, k=1 




<x=1/2 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



T 2 (a + 1) 



r a D?(r a D? R) + k 2 



r 2 (a + l 



■R-u z R=0. 

The first two equations are simple. The third equation 
is a fractional generalization of the Bessel equation [n], 
ll2T ] . It can be solved by fractional series expansion. Since 
Bessel equation is singular at r — 0. We must use such 
ansatz: R = r ap J2m=o c mf am - Substitute it into the 
above equation, we get 



Co 



Cl 



r(q/j + l) 
T(ap — a + 1) 

r(qp + a + l) 
T(ap+1) 



( Tjap_ 



am + 1) 



^ 2 r 2 (a + i) 
i, 2 r 2 ( a + i) 
-^ 2 r 2 («+i) 



= 0, 



\T(ap + am — a + 1 
To have a starting term, cq ^ 0, so 



(30) 
= 0, (31) 



-c m - 2 k 2 = 0. 



( r(ap + l) 
\T(ap-a + l) 



-v 2 T 2 (a + l) 



= 0, 



(32) 
(33) 



FIG. 5: Fractional Bessel function. The fractional index a 
shows up in the exponentials of the fractional Bessel series; a 
small difference from 1 changes the profile largely. For a big 
a, p is small, so the change is suppressed. 



and ci = 0. 

By the recursive relation, a solution is implied 



R P (r) 



= r ap ^{-\) n d n k 2n r a - 2n , 

n=0 



where do = 1 

dn — dn — 1 - 



1 



( r(ap+a-2n+l) 
1 r(ap+a-2n-a+l) 



^ 2 r 2 (a + i) 



(34) 



(35) 



and p satisfies Eq. ([33|) . 

Eq. ([3^|) in one-order calculus (a = 1) is simply p = ±v. 
We meshed in Fig. [4] the surface defined by the equation 
Q33[) . After solving the equation with v = 3 for p when a 
varies, we plotted in Fig. \5\R p (r) belonging to different 
values of a. 
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FIG. 6: Fractional confluent hypergeometric function. 



Thus we get a solution of the above differential equation, 



y( z ) = 



^ {c)% T(ka + iy 



(40) 



Here (a)£ is defined as 
« = 1, « = a, 

r(o + i) 



(o)f a 



r(i) 



r(fca -a + 1) 



r(fca - 2a + 1)/ ' 
fc > 2. (41) 

This can be seen as a fractional generalization of the 
rising factorial 



V. FRACTIONAL HYPERGEOMETRIC 
FUNCTION 

There are other types of special functions in mathemat- 
ical physics. A most famous one is the hypergeometric 
function [l2l - [l4j ]. In this section, we will try to define a 
fractional generalization of the hypergeometric functions. 

Let's first consider the generalization of the confluent 
hypergeometric differential equation: 



z a {D a Yy + (c- z a )D a y - ay = 0. 



(36) 



Here a and c are complex parameters. When a = 1, this 
is the ordinary confluent hypergeometric equation. 
Introducing the fractional Taylor series 



k=0 



(37) 



and substituting, we get the ratio of successive coeffi- 
cients 



Cfc+i • T(ka + a + 1) 
Ck ■ T(ka + 1) 

Cl -r(a + i) 



r(fcq+l) 
r(fca— a+1) 

r(fca+l~ 
r(fca-a+l) 



CO 



a 
c 



(38) 



(39) 



(a)fc = a(a + l)...(a + A; — 1). 



(42) 



And the series (|40|) can be seen as a fractional generaliza- 
tion the confluent hypergeometric function. If a = 1, it is 
exactly the confluent hypergeometric function. Profiles 
of this series (|40l) with different values of a are displayed 
in Fig. [5] 

For the fractional Gauss hypergeometric function, con- 
sider the following series 



y(z) 



mm 



^ [c)% r(ka + iy 



(43) 



which reduces to the Gauss hypergeometric series when 
a = 1. 

The ratio of successive coefficients is 



c k+ i ■ T(ka + a + 1) 
Ck ■ T(ka + 1) 



r(fcq + l) 
r(fea-a+l) 



6 + 



r(fcq+l) 
r(fea-Q+l) 



. I r(fcq+l) \ 
" r r(fca-o+l) I 



or 



(44) 



IYfca + l) T(fca + a + l) 

Cfc+i • c— r - rT + C fe+ l • 



T(ka -a + 1) 



T(ka-a+l) 



, s r(fca + l) T(fca + 1) r(fca + i) 

c k -ab + Ck ■ (a + b)^ ^— + c fc • 



r(fea-a+l) T(fca - a + 1) T(fca - a + 1) 
I 



(45) 



Since 



^V*) = f> ^ (fca + 1) .. z°- fc 



/(z)=^ Cfc z- fc , 



(46) 



fc=i 



' r(fea — a + 1) 



k=0 



(47) 
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FIG. 7: Fractional Gauss hypergeometric function. Gauss 
hypergeometric function is divergent at 1. From the figure we 
can see this is not the case for fractional equivalences with a 
bigger a. 

the equation (|45|) can be translated to a fractional differ- 
ential equation 

ab ■ y(z) + (a + b)z a D a y{z) + z a D a [z a D a y(z)] 

= c-D a y(z) + z a (D a ) 2 y(z). (48) 



When a = 1 , this equation reduces to the ordinary Gauss 
hypergeometric equation. We draw curves of some exam- 
ple functions in Fig. [7] 



VI. SUMMARY 

In this letter, we defined fractional vector calculus in 
the spherical coordinate framework. We discussed Lapla- 
cian equation and heat conduct equation in this kind of 
framework. Some special functions are generalized. 

The geometry induced by fractional operations is non- 
addictive and nonlocal. This kind of geometry is not 
Riemannian; since fractional calculus has found applica- 
tions in many areas of science and engineering, it is the 
effective geometry of such physics processes. It will be 
meaningful to investigate further this kind of geometry. 

For a complete fractional vector calculus, fractional 
vector integral is indispensable. But in this letter we re- 
strained ourselves from this direction. It will be interest- 
ing to discuss fractional vector integral in our framework 
and generalize Green's, Stokes' and Gauss's theorem. 

Special functions show up in different areas of physics 
and engineering. A fractional generalization of these 
functions may find applications in similar situations with 
anomalous tailing behaviors and/or nonlocal properties. 
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